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We present in detail our developed model [Saleh et aL, Phys. Rev. Lett. 107] that governs 
pulse propagation in hollow-core photonic crystal fibers filled by an ionizing gas. By using per- 
turbative methods, we find that the photoionization process induces the opposite phenomenon of 
the well-known Raman self- frequency red-shift of solitons in solid-core glass fibers, as was recently 
experimentally demonstrated [Holzer et a/., Phys. Rev. Lett. 107]. This process is only limited 
by ionization losses, and leads to a constant acceleration of solitons in the time domain with a 
continuous blue-shift in the frequency domain. By applying the Gagnon-Belanger gauge transfor- 
mation, multi-peak 'inverted gravity-like' solitary waves are predicted. We also demonstrate that 
the pulse dynamics shows the ejection of solitons during propagation in such fibers, analogous to 
what happens in conventional solid-core fibers. Moreover, unconventional long-range non-local in- 
teractions between temporally distant solitons, unique of gas plasma systems, are predicted and 
studied. Finally, the effects of higher-order dispersion coefficients and the shock operator on the 
pulse dynamics are investigated, showing that the resonant radiation in the UV [Joly et a/., Phys. 
Rev. Lett. 106] can be improved via plasma formation. 
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INTRODUCTION 

The invention of photonic crystal fibers (PCFs) has led 
to a revolution in the field of nonlinear fiber optics [1^. 
Hollow-core PCFs (HC-PCFs) with the so-called kagome- 
lattice claddings have recently become a superior host for 
the investigation of light-matter interactions between in- 
tense ultrashort optical pulses and gaseous or liquid me- 
dia [2]. These fibers are characterized by a broadband 
transmission range with low group velocity dispersion 
(GVD), and a high confinement of light in the core [3]. 
Many important phenomena have already been observed 
in HC-PCFs. For instance, a drastic reduction in the Ra- 
man threshold has been shown in H2-filled HC-PCFs [4]. 
HC-PCFs have also been used to generate solitary pulses 
by backward stimulated Raman scattering [5 and to ob- 
serve self-similarity in the evolution of transient stimu- 
lated Raman scattering [6 . Important applications such 
as high-harmonic [7 and efficient deep UV generation 
[8] have been recently demonstrated in HC-PCFs filled 
with noble gases, which do not suffer of the limitations 
introduced by the Raman effect. 

The concept of soliton self-frequency blue-shift has 
been introduced and predicted in [9 . In tapered solid- 
core photonic crystal fibers, soliton blue-shift has been 
observed due to the variation of the zero-dispersion wave- 
length (ZDW) along the fiber [10 . In the photoionization 
regime, a limited (of only a few nanometers) ionization- 
induced blue-shift of guided ultrashort pulses has been 
predicted in conventional gas-filled photonic bandgap 
HC-PCFs, which suffer of large group velocity dispersion 
variations near the band edges, thus preventing the clear 
formation of solitary waves [11] |T2] . Very recently, how- 
ever, a sequence of emitted strong blue-shifted pulses has 



been observed in a groundbreaking experiment involving 
Ar-filled kagome-style HC-PCFs, by using few /iJ fs laser 
pulses [13^ . The use of the kagome HC-PCF in these 
experiments was essential due to its unique guiding fea- 
tures. Such fibers possess an unusual broadband guid- 
ance, and a remarkably small group velocity dispersion 
(GVD) (1/^2 1 < 10 fsVcm = 1 psVkm from 400 to 1000 
nm) [HI HH [15J in comparison to the traditional solid-core 
fibers. Moreover, the gas and waveguide contributions 
to the GVD can be balanced in the optical wavelength 
regime, unlike large-bore capillary-based systems, where 
the normal dispersion of the gas always dominates over 
the anomalous dispersion of the waveguide [T5] . 

Prior to our present work on the optical nonlineari- 
ties induced by photoionization, Geissler et. al. modeled 
the photoionization process in terms of the full electric 
field of the pulse [16 . Recently, we have presented a new 
model based on the evolution of the complex pulse en- 
velope [17 , which allowed several interesting predictions 
and it is more suitable for analytical and numerical in- 
vestigation than other previously proposed models. The 
validity of the model has been carefully verified by us- 
ing a numerical model [18 based on the unidirectional 
wave equation [19 . Our model allows us to use well- 
known analytical and numerical methods, developed to 
study pulse propagation in optical fibers, such as per- 
turbation theory, the split-step Fourier method [20] , and 
soliton particle-like analogy [2TJ [22] . By using our new 
model: (i) we derive analytical expressions for the soli- 
ton self-frequency blue-shift in the presence and absence 
of the photoionization threshold of the gas, incidentally 
introducing the completely new concept of floating soli- 
tons; (ii) we predict the possibility of observing two-peak 
'inverted gravity-like' solitons; and finally (iii) we ex- 
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plore the effects of higher-order dispersion coefficients 
and shock term, demonstrating that a plasma-assisted 
blueshift can enhance the efficiency of resonant radiation 
in the UV, possibly improving previous experimental re- 
sults . Our theoretical results closely follow the exper- 
imental measurements in [14 . 

The paper is organized as follows. In Sec. II, the 
governing equations of pulse propagation in HC-PCFs 
filled with an ionizing gas are reviewed. In Sec. Ill, 
the perturbation theory is applied to study the effect of 
photoionization on the soliton amplitude, the temporal 
location of the soliton peak, and the self-frequency shift. 
The possibility of observing two-peak 'inverted gravity- 
like' solitons is discussed in Sec. IV. Section V is dedi- 
cated to study soliton dynamics in plasma. The effects 
of higher-order dispersion coefficients and shock term on 
pulse propagation are examined in Sec. VI. Our conclu- 
sions are provided in Sec. VII. 



GOVERNING EQUATIONS 
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FIG. 1. (Color online), (a) Comparison of the dependence 
of the argon ionization rate on the pulse intensity using the 
full model of Eq. ([T]) and the hnearized model Eq. (b) 
A sketch that represents the attenuation of an initial pulse to 
a floating pulse above the threshold limit due to ionization- 
induced loss. 



Photoionization can take place by either tunneling or 
multiphoton processes. These regimes are characterized 
by the Keldysh parameter pk [ISl [24] . For optical pulses 
with intensities in the range 100 TW/cm^ the Keldysh 
parameter is in the limiting regime between the tunnel- 
ing or multiphoton processes. In this case, the Yudin- 
Ivanov modification of the Perelomov- Popov- Terent'ev 
(PPT) technique is considered to be the appropriate 
model since it accounts for both processes [2F, 26 . Ex- 
perimental measurements [27[ [28 show that tunneling 
models [29l |30] can well describe photoionization pro- 
cesses in noble gases for optical-pulse intensities in the 
range of TW/cm^. This is also confirmed by recent cal- 
culations [18] based on the tunneling Ammosov-Delone- 
Krainov (ADK) model [30 that reproduce simultaneous 
experimental results [14 . Hence, the time-averaged ion- 
ization rate W(/) can be expressed as 



(1) 



W(/)=ci /-i/4exp(-C2/-i/2), 



where c\ and are constants, and / is the laser pulse 
intensity. As shown in Fig. [ija), Eq. ([T]) predicts 
an ionization rate that has an exponential-like behavior 
for pulse intensities above a certain threshold value /th- 
lonization-induced loss that is due to the absorption of 
photons during plasma formation is proportional to the 
ionization rate [IJ. As a consequence, pulses with large 
intensities well above the threshold limit will have their 
intensities strongly driven back to near /th, where the 
ionization-induced loss is drastically reduced. This al- 
lows us to linearize Eq. ([T]) by using the first-order Tay- 
lor series in the proximity of /th, where pulses can survive 
without appreciable attenuation for relatively long time. 
This is our concept of floating pulses introduced in [17] 



and shown in Fig. [TJb). Expanding Eq. ([T]) in its linear 
regime results in 



« a (J - Jth) e (J - Jth) , 



(2) 



where /th and a can be related to ci, C2, and the expan- 
sion point, which is chosen to reproduce the physically 
observed threshold intensity of the considered gas [17]. 
The purpose of the Heaviside function O is to cut the 
ionization rate below Ith- As shown in Fig. [ija), the lin- 
earized model underestimates the ionization rate thence 
the ionization loss, in comparison to the full model. This 
yields surprisingly to a similar qualitative behavior be- 
tween the two models even for / > /th, since the ioniza- 
tion rate and the ionization loss are the key factors in the 
photoionization process. 

By using Eq. one can show that propagation of 
light in a HC-PCF filled with an ionized Raman-active 
gas can be modeled by the following coupled equations: 



id, + D{idt)+jKR{t)(E)\'^{t)\^ 



2fcoc2 



= 0, 



(3) 

dtn, = [<?Meff] [nx - ne] - 9 - I*]?,) , 

(4) 

where ^(^,t) is the electric field envelope, z is the lon- 
gitudinal coordinate along the fiber, t is the time in 
a reference frame moves with the pulse group velocity, 
D{idt) = Xlm>2 Pmi'^dt)^ /ml is the full dispersion oper- 
ator, /3m is the mth order dispersion coefficient calculated 
at an arbitrary reference frequency cjq, 7k is the Kerr 
nonlinear coefficient of the gas, R{t) = (1 — p) S{t)-\-ph{t) 
is the normalized Kerr and Raman response function of 
the gas, S{t) is the Dirac delta function, p is the rela- 
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tive strength of the non- instantaneous Raman nonhnear- 
ity, h{t) is the causal Raman response function of the 
gas [TTJ |20] , the symbol (g) denotes the time convolution 
{A^B = jA{t- t')B(t')dt' = jB{t- t')A(t')dt') , c is 
the speed of light, /cq = cjq/c, uoq is the input pulse cen- 
tral frequency, uj^ = [e^ne/(eome)]^^^ is the plasma fre- 
quency associated to an electron density ne(t), e and ttIq 
are the electron charge and mass respectively, and eo is 
the vacuum permittivity, a = ai -\- a2 is the total loss 
coefficient, ai is the fiber loss, a2 = dtfiQ is the 

ionization-induced loss term, Aeff is the effective mode 
area, Uj is the ionization energy of the gas, = /^eff, 
l^l^j^ = /th^eff, and nx is the total number density of 
ionizable atoms in the fiber, associated to the maximum 
plasma frequency ujt = [e^nT/(eo^e)]^^^- The recombi- 
nation process is safely neglected in these coupled equa- 
tions since the pulse duration (of the order of tens of fs) 
is always shorter than the recombination time [24l [32] . 
However, the recombination effects as well as the rear- 
rangement of the electronic plasma due to ponderomo- 
tive forces should be included in the case of long-pulse 
propagation, making the problem considerably difficult 
[331 [34 . 

Introducing the following rescalings and redefinitions: 

^ = Z/Zo, T = t/to, ^0 = [7K^0]"^/^ = ^/^O, 

IV^thl = r(r) = R{t)to, (j) = \kQZQ[uj^/ujQf , 

= ^/co^^o [cjt/^o]^, cf = ato/[^eff7K^o], and k = 
[// a eo me cc;o/[/co e^], where zq = to/l/^2('^o)| is the 
second-order dispersion length at the reference frequency 
cjo, and to is the input pulse duration [20]. Hence, the 
two coupled equations for fioating pulses can be replaced 

by 



id^ + D{idr)+r{T)®\tp{T)\'^ -(l) + ia V = 0, (5) 
drd^ = a{<i>^ - <i>) - e - |Vlt\) , (6) 



(ODEs) are obtained: 



where a = k {4>t - <j>) [l - \i^\lJ\i^?] 6 - IV'I 
and the fiber losses are neglected. 
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PERTURBATION THEORY FOR SOLITONS 

The effect of the Raman and ionization perturba- 
tions on the soliton dynamics can be studied using 
the perturbation theory described in standard text- 
books (e.g. [20]). Neglecting higher-order disper- 
sion coefficients, i.e. Pm>2 = 0, the solution of 
a perturbed nonlinear Schrodinger equation, id^ip + 
h d^ilj + IV^PV^ = ^e(?/^) is assumed to be ?/^(^,r) = 



A(e)sech[A(C)(7 



where A{^) is the 



soliton amplitude, rp(^) is the temporal location of the 
soliton peak, is self- frequency shift, and e{ip) is 

the perturbation function. Substituting this Ansatz in 
Eqs. (5]6), the following ordinary differential equations 



= -Im|y" e(^/;)tanh[A(r-rp)] ^/;*dr| , (7) 
^Re||_^"e(V')(T-rp)V*5T|, (8) 

(9) 



d5 

A 

OA ( /'"'"°° "1 
^ = Re 

where Re and Im stand for the real and imaginary parts. 

In order to extract useful analytical information from 
Eqs. ( [5j6l ), we will start with the simple case where 
we assume that the ionization loss and the thresh- 
old are negligible. Ignoring ionization loss is unrealis- 
tic unless the pulse under consideration has its max- 
imum amplitude just above the threshold, in which 
case it does not feel strong plasma-induced decay. 
Hence, Eq. (|6| can be solved analytically, 0(r) = 



(/)T 1 1 — exp 
dition (j){—oo) 



~^ S-oo \i^{^')\^dT' |, with an initial con- 
= 0, corresponding to the absence of any 
plasma before the pulse arrives. For a small ionization 
cross-section, ^(r) r] |?/^(r') pdr', where r] = 
Moreover, in the long-pulse limit \iI){t — t')\^ |?/^(r)p — 
r^dr\il^{r)\'^ [20 . This allows the two coupled equations 
to be reduced to a single partial integro-differential equa- 
tion: 



(10) 

where tr = r' r{r') dr' . This equation shows clearly 
that the effect of ionization is almost exactly opposite to 
that of the Raman effect: the fourth term in Eq. (10) 



involves a derivative of the field intensity, while the fifth 
term involves an integral on the same quantity. Hence, 
one can conjecture that the ionization perturbation will 
lead to a soliton self-frequency blue-shift, instead of a 
redshift. However, this is not an exact analogy because 
of thermodynamical reasons. In fact, the Raman ef- 
fect acts on a soliton by constantly decreasing its cen- 
tral frequency, leaving its pulse shape in the time do- 
main and in the frequency domain undeformed - in other 
words the soliton continuously converts its coherent en- 
ergy into incoherent optical phonons in the medium, and 
therefore the number of photons contained in the soli- 
ton does not change, and the overall entropy increases 
according to the second law of thermodynamics. How- 
ever, in the photoionization-induced soliton blueshift de- 
scribed above, the soliton coherently receives energy from 
the medium (the plasma in this case). For this to be 
physically possible without breaking the second law of 
thermodynamics, the number of photons in the soliton 
must necessarily decrease - thus the unavoidable plasma- 
induced losses cannot be eliminated on principle. The 
stronger the initial pulse intensity, the stronger the soli- 
ton blueshift, and the larger the plasma-induced losses. 
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so that Eq. (10), which assumes smah losses, would lose 



its validity. However, in the following we shall prove nu- 
merically that when solitons decrease their amplitudes to 
just above the ionization threshold, their losses are very 
limited thence they can propagate for relatively long dis- 
tances. This led us to naturally introduce the concept of 
floating solitons^ i.e. those solitary pulses with maximum 
amplitude just a little above the threshold — the regime 
where ionization losses can be neglected and Eq. ( [1Q| ) 
maintains its validity [17 . 

Solving Eq. ([To]) with a perturbation function e(?/^) = 

-i^[mdrm^^r]J^^mr')\^dr'], results in A(0 = 

A{0) = Ao, 5{0 = ^Raman(0 + ^ion(0 = '9 ^p(0 = 

and g = g^ed + ^biue, where ^red = +(8/15)rRA^ 
and ^biue = "(2/3)7^^0 P^- Note that g can be posi- 
tive, negative or even zero, depending on the value of 77, 
tr and Aq. The precise rate for the self- frequency blue- 
shift can be obtained using the exact formula of ^(r). 



= cr-Mo Vt [(1 - crAo) - (1 + aAo) exp(-2crAo) 



which tends to ^biue foi" small values of <j, but start to 
differ considerably from it for Aq > <j~^. 

A further step can be applied to take into account the 
effect of the threshold intensity on the blue-shift rate. 
For floating pulses, the ionization loss is not large and 
can be neglected to a good approximation. For such 
pulses, only a small portion of energy above the thresh- 
old intensity contributes to the electron density build-up. 
In this case, 0(r) o:^ r] J^rf^ [|?/;(r')p — IV^I^h] dr^ where 
T = Aq ^sech~^ [|?/^|th/Ao]- This formulation embeds the 
Heaviside function introduced in Eq. (|6|. Replacing the 
fifth term in Eq. ([To| with the new definition of ^(r) and 
solving Eq. ([7|), results in 



-riAl 



tanh'' e 



1^1 



^0 



{esech^e 



tanh6') 



(11) 

with = AqT. It can be seen that this expression tends 
to -{2/3)r]Al for smah |V;|th. 

The threshold intensity effect on the blue-shift rate can 



also be alternatively estimated by writing Eq. ( 10 ) as 



J —00 

(12) 

where fj = r] fi and 
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FIG. 2. (Color online). The spatial dependence of the soli- 
ton amplitude and frequency shift during a photoionization 
process for different initial pulse amplitudes. 



Moreover, the effect of the photoionization loss on 
the soliton amplitude and frequency shift can be stud- 
ied numerically by solving the above ODEs with e(?/^) = 



where T is de- 



_T mr'r-\i'\i^)dT' 
termined via its previous definition using the initial pulse 
amplitude, and the variation of the position T due to 
losses is neglected, 



dA 



-2k (^T (Atanht? - Itp] 



th - 



rjA'^ 



tanh"^ i9 



th 



(^sech^i?- tanh^) 



(14) 

= AT. These equations can be solved numerically to 
determine the spatial dependence of the soliton ampli- 
tude and frequency shift, as shown in Fig. [2] Pulses 
with initial large intensities (Aq > IV^Ith) have a boosted 
self- frequency blue-shift. However, the ionization loss 
suppresses the soliton intensity after a short propagation 
distance to the float ing-soliton regime, where the soli- 
ton can propagate for a long propagation distance with a 
limited blue-shift and negligible loss. The maximum fre- 
quency shift is achieved when the soliton intensity goes 
below the photoionization threshold. 



/5t 0^(^01' -l^lth] dr' 



(13) 



The factor fi represents the ratio between the pulse en- 
ergy contributing to plasma formation and the total en- 
ergy of the pulse, see Fig. [ijb). Hence, it takes into 
account the overestimation of the plasma density by ne- 
glecting the threshold intensity. As a consequence, one 
can simply prove that ^biue = ~{2/3)vAq. 



INVERTED GRAVITY-LIKE BOUND STATES 

The previous conclusions on the strong analogies be- 
tween photoionization shift and Raman shift of solitons 
allow us to construct a novel kind of solitary wave that 
can be supported by both the Kerr and the ionization 
nonlinearities. This solitary waves will be analogous to 
the Raman bound solitons discovered experimentally in 
^35j and explained theoretically in [36H38] . In the non- 
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inertial reference frame of an accelerated soliton (with 
a new time coordinate ( = r — ^^^/2), and by us- 
ing the so-called Gagnon-Belanger gauge transformation 
HC,r) = fiOexv[i{q-g^e/3 + 9T)^] [39 , Eq. (|l2) 
can be written as an integro-differential equation, 



dif 



UiC)f = qf, 



(15) 



where ionization loss and higher-order dispersion are ne- 
glected, UiO =9C~ l/P + rud^lfl^ + i) j^^ \f\^dC is 
a potential [40l |4T], in which the first term is gravity- 
like, q = Aq/2 — fjAo is the soliton wavenumber, and Aq 
is the amplitude of the one-peak solitary solution of Eq. 



(15). In the case when only Raman effect was present. 



a special kind of nonlinear metastable bound states have 
been found as a general solution for Eq. (15) [4QH42]. 



Such bound states are multi-peak stationary states (in 
the non-inertial reference frame moving with acceleration 
9red)^ which are due to the Kerr effect complemented by 
the Raman nonlinearity [36ti38] . 

In the presence of a Raman- inactive gas (such as Ar- 
gon) inside the HC-PCF, we have seen that the ioniza- 
tion process leads to a soliton acceleration in the time 
domain, and to a linear shift of its frequency toward the 
blue with a rate ^biue- In this case, solitons will feel an 



' ant i- gravity' field, and the solution of Eq. (15) can be a 



multi-peak stationary state with a negative slope oppo- 
site to the Raman case. An example of a two-peak bound 
state found numerically by using the shooting method is 
depicted in Fig. [sja). Multi-peak solitary solutions anal- 
ogous to those found in |37j can also be obtained. The 
propagation of the two-peak bound state found in Fig. 
Ija) in a long HC-PCF is shown in Fig. ^h). 

The characteristics of a two-peak solitary solution — 
the temporal separation ("q and the amplitude imbalance 
R between the two peaks — can be analytically deter- 
mined. Consider a general two-peak solitary solution, 

/(C) = Aisech [A, (C - Ci)] + ^2sech [A^ (C - C2)] , (16) 

where A2 = RAi and Co = C2 — Ci- principle, a set of 
three algebraic equations is needed to determine Ai, Co 



and R. By substituting the above solution in Eq. (15), 
the first equation is obtained. 



^biueCo = Al{l - R^)/2 + 7)Ai(l + R), 



(17) 



where fj is assumed to be the same for the two solitons. 
The last term on the right-hand side is due to the inte- 
gration constant in the nonlinear potential U{Q, and it 
is not present in the analogous expression for the pure 
Raman-bound states [i.e. — ^redCo = Al{l — R'^)/2] since 
the Raman effect is modeled in the equations by a deriva- 



tive. Multiply Eq. ( 15 ) by /* and integrate its both sides, 
results in 



-/*a|/-/7(c)|/p 



dc 



\f\'dC (18) 
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FIG. 3. (Color online), (a) A two-peak soliton stationary 
solution for a Raman- inactive gas with f) = 0.02 and q = 
0.3. The dashed-dotted red, dashed blue, and solid green 
curves represent the nonlinear potential U, the pulse intensity 
and the ionization field 0, respectively, (b) Temporal 
evolution of the two-peak soliton in a gas-filled HC-PCF fiber. 



where the superscript * denotes the complex conjugate. 
A second equation can be attained by substituting Eq. 
^ into Eq. (ITs]), 



Al{l^R^)/2-f]Ai{1^2R^R^)-q{l^R) ^ 0. (19) 
The set of three algebraic equations is completed by Eq. 



(15) evaluated at C = 0. Hence, the variables Ai, Co and 
R can be determined for a certain value of q. 

Finally in Raman-active gases, solitons may feel ei- 
ther a 'gravity- like' or an 'inverted- gravity-like' poten- 
tial based on the sign of the total gravity acceleration 
g = ^red + ^biue- Hcucc, the two-pcak solitary solution 
can have either positive or negative slopes as depicted 
in the panels (a,b) of Fig. [4j respectively. Two-peak 
solitons with positive (negative) slope are obtained when 
Raman (photoionization) process is dominant. 



SOLITON DYNAMICS IN PLASMA 

In this section, we will elaborate our results in [17 con- 
cerning soliton dynamics in an ionized gas. Equations 
( 5j6 ) can be numerically solved by using the split-step 
Fourier method [20] in order to fully describe the prop- 
agation of pulses in the presence of plasma. In this sec- 
tion, higher-order dispersion coefficients are assumed to 
be negligible, i.e., /3m>2 = 0. The temporal and spectral 
evolution of a higher-order sech-pulse, with an initial in- 
tensity less than the threshold value, are depicted in the 
panels (a,b) of Fig. [sj respectively. Panel (c) shows the 
variation of the ionization fraction along the fiber. An 
optical pulse pumped in the deep anomalous-dispersion 
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FIG. 4. (Color online). A two-peak soliton stationary so- 
lution in a Raman-active gas with fj = 0.02 and q = 0.3. 
(a) Gravity-like bound states with tr = 0.1. (b) Inverted 
gravity- like bound states with tr = 0.01. The dashed-dotted 
red, dashed blue, and solid green curves represent the non- 
linear potential [/, the pulse intensity and the ionization 
field (/), respectively. 



regime of the fiber undergoes self-compression. When 
the pulse intensity exceeds the threshold value, a certain 
amount of plasma is generated due to gas ionization, and 
a fundamental soliton is ejected from the input pulse. 
The soliton central frequency continues to shift towards 
the blue due to the energy received from the generated 
plasma. However, due to the concurrent ionization loss, 
the soliton intensity gradually decreases until it reaches 
the threshold, thence the pulse blue-shift is ceased. A 
second ionization event accompanied by a second-soliton 
emission can take place by further self-compression of the 
input pulse based on its initial intensity. 

A clear representation for the pulse dynamics in plasma 
is shown in Fig. [6j where the temporal profile of the 
pulse intensity is plotted at selected positions inside 
the fiber. The simulation parameters are similar to Fig. 
[5] Initially, |?/^(r)p is insufficient for plasma ionization. 
Due to soliton-breathing, the pulse passes through a self- 
compression stage which strongly enhances its maximum 
intensity [43^. The amount of optical energy above the 
threshold intensity l^l^h contributes to plasma formation 
that emits a blue-shifted soliton. However due to the 
ionizat ion-induced loss, the soliton amplitude is atten- 
uated to the regime where |'0(t)P ^ IV^Ith- Such pulses 
are dubbed floating solitons^ since they can propagate for 
considerably long distance with minimal attenuation and 
limited blue-shift. Based on the initial input-pulse inten- 
sity, other solitons can also be emitted due to further 
self-compression. At the end, a train of floating solitons 
are generated. Indeed, Fig. |6]clearly shows the formation 
of a series of floating solitons, that exist and propagate 
for relatively long distances with their maximum ampli- 
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FIG. 5. (Color online). Temporal (a) and spectral (b) evo- 
lution of an energetic pulse propagates in an Ar-filled HC- 
PCF. The temporal profile of the input pulse is A^sechr, with 
N = 8. The gas pressure is 5 bar. The panels show the 
ejection of two solitons that continue blue-shifting until the 
ionization loss suppresses their intensities below the threshold 
value. Contour plots in this paper are given in a logarithmic 
scale, (c) Spatial dependence of the ionization fraction along 
the fiber. 



tude just above the ionization threshold. These objects 
are completely unknown in conventional solid-core fiber 
optics. The fact that these results have been confirmed 
in concurrent experiments [M] is a very convincing proof 
of the validity of our master equations (5)6). 

We have also found an interesting non-local interac- 
tion between successive solitons due to the non- vanishing 
electron density tail, when their temporal separation is 
shorter than the recombination time. Due to this inter- 
action, a leading soliton can slow down the acceleration 
of a trailing soliton by an exponential factor. In the fre- 
quency domain, the leading soliton suppresses exponen- 
tially the blue-shift of the trailing soliton. The reason is 
that the ionization field ^(r), created by the first soli- 
ton, decays at a relatively slow rate. This establishes a 
unique 'non-local' interaction between this soliton and 
other temporally distant solitons. 

These unprecedented dynamics are featured in Fig. [7| 
that show the temporal and spectral dependence on the 
soliton order N assuming that the input pulse is sechr. 
A clustering between two or more distant solitons in both 
temporal and spectral domains is also observed at some 
'magic' input energy as a result of the interplay between 
the non-local interaction, ionization loss and ionization 
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FIG. 6. (Color online). Intensity profile of a pulse in the time 
domain at different positions, ^, inside an Ar-filled HC-PCF. 
The dashed red line represents the threshold intensity. The 
simulation parameters are similar to Fig. [5] Each panel is 
titled by its main feature. Insets are enclosed in panels for 
better view and more details. 



FIG. 7. (Color online). Dependence of the temporal (a) and 
spectral (b) outputs of an energetic pulse Nsechr on the 
soliton order N. The fiber is an Ar-filled HC-PCF with length 
^ = 1/4. The gas pressure is 5 bar. Temporal and spectral 
clustering occur at N = 9.2 due to the long-range non-local 
soliton interactions, (c) The output energy versus the input 
energy. 



threshold. The scenario is as follows: As long as the 
intensity of the first-emitted soliton, /i, is above the 
threshold intensity /th, it prevents the ejection of a sec- 
ond soliton due to the strong effect of the non-local in- 
teraction. As soon as, Ii approaches /th by virtue of the 
ionization loss, the first soliton becomes a floating soli- 
ton. Hence, the blue-shift and acceleration of this soli- 
ton are reduced significantly. Simultaneously, the second 
soliton is emitted and it can recover its expected accel- 
eration and blue-shift due to the nearly disappearance of 
the first-soliton non-local force. This allows the second 
soliton to catch and cluster with the first soliton. In ad- 
dition, the spectrum of the two solitons start to overlap 
and form spectral clustering. Similarly, the intensity of 
the second soliton approaches /th due to the ionization 
loss, allowing a third soliton to cluster with the other two 
solitons. When the first two solitons are very close to 
each other, they push back the third soliton due to their 
combined non-local force. In fact, the dynamics after the 
clustering becomes too complicated to be interpreted in 
simple terms. Figure [tJ^c) shows the output energy versus 
the input energy. The linear dependence of the output 
energy is slightly broken at the points, which correspond 
to the ejection of a new soliton. Our results show an ex- 
cellent qualitative agreement with other non-analytical 
numerical techniques [18 . 

XFROG spectrograms for pulses with initial temporal 
profile A^sechr and different initial intensities are de- 
picted in the panels of Fig. [sj where (a) represents the 
reference pulse; (b) and (c) shows the emission of the 
first and second solitons, respectively; and (d) depicts 
the temporal and spectral clustering of the first two soli- 
tons and the emission of a third soliton. 
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FIG. 8. (Color online). XFROG spectrograms for pulses with 
selected soliton order N — in an increasing order. The simu- 
lation parameters are similar to Fig. [t] (a) = 3, ^ = 0. (b) 
AT = 5, ^ = 1/4. (c) AT = 8, ^ = 1/4. (d) N = 9.2, ^ = 1/4. 
Each panel is titled by its main feature. White arrows show 
the movement of the solitons. 



GENERATION OF DISPERSIVE WAVES 

The contribution of the higher-order dispersion coeffi- 
cients Prn>2 to the pulse dynamics starts to play a signif- 
icant role as pulse central wavelength moves towards the 
zero dispersion wavelength (ZDW), where /32(a;o) ^ 0. 
In fact, higher-order dispersion coefficients may lead to 
a phase matching condition between two different waves 
— an optical pulse, with a central wavelength lies in the 
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FIG. 9. (Color online). Spectral evolution of an energetic 
pulse A/'sechr, with A/" = 10, under the influence of higher- 
order dispersion coefficients with: (a) switching off the ion- 
ization process, (b) switching on the ionization process. The 
optical pulse central wavelength is 0.8 /im and the gas pressure 
is 2 bar. The dotted line represents the ZDW. 



FIG. 10. (Color online). Spectral evolution of an energetic 
pulse under the influence of higher-order dispersion coeffi- 
cients and the shock term with: (a) switching off the ion- 
ization process, (b) switching on the ionization process. The 
simulation parameters are similar to Fig. [9] except that the 
shock operator Tghock is included. The dotted line represents 
the ZDW. 



anomalous-dispersion regime and is close to the ZDW, 
and a dispersive resonance wave in the normal-dispersion 
regime [HI 05] . Recently, this fact was implemented to 
achieve a coherent deep-UV laser source by using an Ar- 
filled HC-PCF [8]. Joly et al have obtained 8 % con- 
version efficiency from IR to deep-UV, where the output 
can be tuned easily via the pulse energy and gas pres- 
sure. The pulse power level was kept below the ionization 
threshold to avoid any ionization-induced loss. However, 
we find an enhancement in the dispersive- wave radiation 
via the photoionization process by more than an order of 
magnitude due to the ionization-induced blue-shift. As 
the optical pulse moves towards the blue or the ZDW, 
new dispersive waves are generated. The conversion to 
the normal-dispersion regime is ceased when the opti- 
cal intensity reaches the threshold value by the ioniza- 
tion loss. Unlike the Raman process, where the pulse 
central wavelength is shifted away from the ZDW dur- 
ing pulse propagation results in arresting the optical fre- 
quency conversion. The effect of the photoionization on 
enhancing the generation of the dispersive waves is shown 
in Fig. [9J The top (bottom) panel represents the spec- 
tral evolution of a higher-order sech-pulse in a gas-filled 
HC-PCF, where the ionization is switched off (on). 

The dispersion of Kerr nonlinearity is usually asso- 
ciated with effects such as self-steepening and optical 
shock formation [20^. Due to the absence of Raman ef- 
fect in kagome-style HC-PCFs filled with noble gases, 
and the weakness of the group velocity dispersion (GVD) 
in such fibers, the role of the Kerr nonlinearity disper- 
sion assumes an unusual importance. This effect can be 
studied by using Eqs. (3j4), where 7k is replaced by 
7k (1 -^iTshockdt), and Tghock = l/^o is the shock time 



[20^. The importance of the shock term is depicted in 
Fig. [lOj where the simulation parameters are similar to 
Fig. [9] except that Tghock is included. Involving the shock 
term in the pulse dynamics increases the conversion effi- 
ciency to the normal-dispersion regime from 3 x 10~^% 
to 5 X 10~^% at the end of the fiber due to the spectral 
asymmetry of the pulse. 



CONCLUSIONS 

We have presented a detailed model based on the evo- 
lution of the pulse complex envelope to study pulse prop- 
agation in gas-filled HC-PCFs under the influence of 
the nonlinear photoionization-induced effects. By ap- 
plying perturbation theory, we show that the photoion- 
ization process represents the exact counterpart of the 
Raman self-frequency red-shift of solitons when their in- 
tensities are slightly above the threshold intensity. Ex- 
pressions of the soliton self-frequency blue-shift are de- 
rived in the presence and absence of the photoioniza- 
tion threshold. Moreover, the influence of the ioniza- 
tion loss on the soliton amplitude and the frequency 
shift is studied. Using the Gagnon-Belanger gauge trans- 
formation, stationary negative-slope two-peak 'inverted- 
gravity-like' solitary solutions are obtained for pulses 
propagate in HC-PCFs filled by Raman-inactive gases. 
However, positive- or negative-slope two-peak solitary 
solutions can also be attained in the presence of a tun- 
able Raman-active gas. The pulse dynamics, obtained 
by using the split-step Fourier method, shows the soli- 
ton emission, break-up and blue-shift. Furthermore, we 
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find unconventional long-range non-local interactions be- 
tween successive solitons due to the non- vanishing elec- 
tron density tail. The interplay between this unprece- 
dented interaction, the ionization-induced loss and the 
ionization threshold yields a spectral and temporal clus- 
tering between distant solitons. Finally, the effects of 
higher-order dispersion coefficients and shock term are 
exploited in generating a strong dispersive-wave radia- 
tion in the normal-dispersion regime of the fiber, assisted 
by the absence of the Raman effect in noble gases. 

We would like to thank John Travers, Philip Holzer, 
Wonkeun Chang, Nicolas Joly and Philip St. J. Russell for 
useful discussions. This research is funded by the Ger- 
man Max Planck Society for the Advancement of Science 
(MPG). 
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